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We derive a formula for the quantum corrections to the electrical current for a metal out of 
equilibrium. In the limit of linear current- voltage characteristics our formula reproduces the well 
known Altshuler-Aronov correction to the conductivity of a disordered metal. The current formula is 
obtained by a direct diagrammatic approach, and is shown to agree with what is obtained within the 
Keldysh formulation of the non-linear sigma model. As an application we calculate the current of a 
mesoscopic wire. We find a current-voltage characteristics that scales with eV/kT, and calculate the 
different scaling curves for a wire in the hot-electron regime and in the regime of full non-equilibrium. 



I. INTRODUCTION 

Quantum interference effects in disordered metals have 
been the subject of intensive investigation for over twenty 
years. For general reviews see The interference of 

the scattered electrical waves in the presence of a random 
potential leads to corrections to the semi-classical 

formula of the electrical conductivity known from the 
Drude-Boltzmann theory. The physical implications of 
these quantum corrections have been extensively dis- 
cussed in the literature mostly for the equilibrium prop- 
erties, for which experimental data were available, al- 
though a number of non-linear electric field effects have 
been predicted in the past (2|j^-^ . 

In contrast, non-equilibrium electrical transport has 
received considerable attention in the field of mesoscopic 
physics. Examples are the transport in quantum dots 
pO| , or the shot noise in mesoscopic conductors ||T^. In 
this situation however the majority of the studied phe- 
nomena did not involve interference effects. 

Our interest in the non-equilibrium properties of in- 
terference phenomena originated by the suggestion fl^ j 
that the non-equilibrium electric noise could be the ori- 
gin for the low temperature saturation of the weak lo- 
calization dephasing time observed in disordered films 
and wires jl^. Mohanty et al. pointed out that in 
the samples with the strongest dephasing rate also the 
interaction correction (Altshuler-Aronov) to the conduc- 
tivity saturates at low temperature. This suggests that 
also the Altshuler-Aronov correction should be affected 
by the non-equilibrium noise. Similar speculations con- 
cerning an electric field effect on the Altshuler-Aronov 
correction have also been made earlier by different au- 
thors in different contexts [ p^Hl^ . 

This motivated us to study the interaction correction 
in the presence of an external field. In analogy to the 
analysis performed in the literature for weak localization 
1^, we calculated the interaction correction to the cur- 
rent in the presence of a time dependent vector poten- 
tial, assuming a thermal distribution function p^ . In- 
deed we verified the above-mentioned speculations since 
we found an electric field effect. However, contrary to 



what is known in the case of weak localization, where the 
strongest effect occurs when the period of the AC field is 
of the order of the dephasing time, we found a suppres- 
sion of the interaction correction even by static electric 
fields. These findings have however raised the issues of 
the possibility of experimentally observing the effect and 
of its physical interpretation. Both these questions will 
be addressed in this paper. 

Let us comment first the problem concerning the inter- 
pretation of the effect. In Ref. ||l^ we demonstrated that 
the non-linear field effect can be understood in terms of 
dephasing by calculating the phase shifts of the relevant 
classical paths in the presence of a time dependent vec- 
tor potential. On the other hand, a static electric field 
can also be described in terms of a static scalar poten- 
tial, and it is clear that a static scalar potential does not 
lead to dephasing. Therefore we think it is of interest to 
present a version of our theory in the scalar gauge. In 
this paper we will explicitly show the gauge invariance of 
our previous results. 

The second problem concerns the scale of the effect and 
its experimental observability. We found in Rcf. jl^ that 
the temperature-dependent Altshuler-Aronov correction 
to the conductivity saturates when the voltage drop on 
the thermal length {Lt ~ y/D/T) is comparable to the 
temperature cELt ~ kT. From this condition one can 
estimate the strength of the microwave field that is nec- 
essary to explain the saturation of the Altshuler-Aronov 
correction observed in the experimental data of Ref. . 
In so doing one arrives at a microwave field value that is 
more than an order of magnitude larger than the opti- 
mistic estimate given in Rcf. to explain the saturation 
of the weak localization time. Furthermore it is rather 
unlikely that the condition cELt ~ kT can be reached at 
low temperature, since strong heating is already assumed 
to set in when the voltage drop over the electron-phonon 
length is of order of the temperature. Strong heating on 
the other hand is not observed in Ref. . Similar prob- 
lems arise also in the attempt to explain the experimental 
data of Refs. [|5|jl6). 

Despite the above mentioned problems, we think we 
cannot rule out non-equilibrium noise as a reason for 
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the observed saturations. In fact, even in the absence 
of strong heating, the distribution function may devi- 
ate from the equihbrium form and affect the interaction 
correction to the conductivity and possibly lead to sat- 
uration at considerably weaker electric fields. A theory 
which is valid even out of equilibrium will be developed 
in this paper. 

In contrast to our previous work we will avoid 
to guess the distribution function which could be rele- 
vant for the experiments of Refs. [ p^p4| . Instead we will 
calculate the interaction correction in a more controlled 
situation. Nowadays it is, indeed, possible to create non- 
equilibrium in a controlled way by, for instance, attaching 
a short mesoscopic wire to large metallic reservoirs (see 
e.g. 10). In the absence of inelastic scattering processes 
the distribution function in the wire is a linear super- 
position of the distribution functions of the leads. The 
interaction correction in such a situation has also been 
considered in the recent paper by Gutmann and Gefen 
|pO| . We verify their result that the I —V characteristics 
scales as eV/kT. Going beyond, we calculate the I — V 
characteristics explicitly and we will compare quantita- 
tively the wire in non-equilibrium with the wire in the 
hot electron regime. In addition we will also discuss the 
interaction correction in the spin triplet channels. 

Our paper is organized as follows. In the next section 
we introduce the basic quantities and recall the main re- 
sults of the Drude-Boltzmann theory within the Keldysh 
formalism. In section III we consider the quantum cor- 
rections to the conductivity within the Keldysh diagram- 
matic approach. We derive, in particular, an expression 
for the current in the presence of an external electric field. 
In section IV we discuss the gauge invariance of the the- 
ory, while in section V we present a specific application: 
a mesoscopic wire. Finally in section VI we give our 
conclusions. In the appendices we outline how to obtain 
the same results using the Keldysh formulation of the 
non-linear sigma model and we extend the calculations 
in order to include also the spin effects. 



G^ix, x') = -ie(t - t') {{<i'{x)<i'\x') + -^\x')^{x))) (2) 
G^ix, x') = +ie(t' - t) {{■^{x)-^^{x') + -^^{x')^ix))) (3) 
G^(x, x') = -i ((*(x)*^(2;')-*^(a;')*(a;))) , (4) 

where \I' and are fermion operators and x — (x, t). In 
equilibrium the Keldysh component of the Green func- 
tion is expressed in terms of the retarded and advanced 
components by Gf = [1 - 2/(e)](Gf - G^), where /(e) 
is the Fermi function. The Keldysh component out of 
equilibrium will be discussed later. 

The Green function solves the differential equation 



d 1 

1— + —(V + ieA)2 + e(j) + fi] G(x, t; x', t') 



dt 2m 



diidxiE(x, t; xi , ii)G(xi , ti; x', t') 

= S{^-^')S{t-t'), (5) 



where and (j), A are the scalar and vector potential. Since 
the self-energy E has the same triangular matrix struc- 
ture as the Green function, one can invert the inverse 
Green function G^^ and finds for the Keldysh compo- 
nent the relation 



G^ = G^E^G^. 



(6) 



For a graphical representation see Fig.^ 
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FIG. 1. Graphical representation of the Green function; 
the shaded box in represents the Keldysh component of 
the self-energy, i.e., basically the distribution function. 

We then introduce the ^-integrated (quasi-classical) 
Green function 



II. BASIC DEFINITIONS AND THE 
DRUDE-BOLTZMANN THEORY 

In this section we will recall some basic relations of 
the quasi-classical approximation in its non-equilibrium 
(Keldysh) formulation [|T| . Our notation will mainly fol- 
low Ref. j2^. We will write down the equation of mo- 
tion for the Green functions in the presence of impurity 
scattering in the case when quantum interference is com- 
pletely neglected. The Green functions have the matrix 
structure 



G - I Q 



(1) 



5„,(P,R)-z / d^clre-'P '-G(R+|,i;R-|,i' 



(7) 



with 



where ^ = p^/2m — /i and p is a unit vector along the 
momentum. The Green function in the energy domain is 



5„'(p,x)= / dMi'e-*--'*'5«'(p,x). (8) 



We will keep the notation of small g for the ^-integrated 
Green functions and capital G for the not integrated 
Green functions all over this paper. When approximating 
the density of states as an energy independent constant, 
the ^-integration is related to an integration over the mo- 
mentum p according to 
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(27r) 



(9) 



We will now recall some relations that are specific for 
impurity scattering. By treating the impurity scatter- 
ing within the self-consistent Born approximation and as- 
suming a Gaussian, 5-correlated impurity potential with 



(L/(x)L/(x')> 
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-(5(x-x'), 



(10) 



2ttNot 

the electron self-energy is local in space and is given by 
1 



where the diffusion constant is I? = VpT/3. Notice that 
this equation is solved by any function g^, (p,x) which 
is independent of position x and which depends on time 
differences {t — t') only. This reflects the fact that any 
distribution function is allowed for noninteracting elec- 
trons. 

The charge density and current density are related to 
the Keldysh component of the Green function, 



E""P(x,t;x',t': 



G'(x,t;x,t')(5(x-x') (11) p(x,<)=ieG^(x,t;x,t) 
(12) 
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s;7(x)j(x-x' 



(19) 



j(x, t) = — [Vx- Vx' +2ieA(x, t)] G^(x, i; x', t)|x'=x. (20) 
Zm 



Notice that this equation has to be solved self- 
consistently for all the components of the Green func- 
tion. Using the above definition, one observes that the 
impurity self-energy is related to the s-wave part of the 
quasi-classical function. 



1 / dp 
2t / 47r- 



(13) 



The distribution function out of equilibrium is found 
by solving the appropriate kinetic equation. Here we de- 
rive the kinetic equation for g^, from Eqs.(^) and @). 
For simplicity we neglect external fields for the time be- 
ing. Under these conditions the retarded and advanced 
Green functions are 
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e-e±i/2r- 



(14) 



Near the Fermi energy (e, e' <^ ep) and for small mo- 
menta {q <C pf) one finds 

ff^-(p,q) = ^/ ^G^(e,p + q/2) 

xI]^,(q)G^(e',p-q/2) (15) 

-gi^,(p,q). (16) 
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r e — e' + i/r — vpp • q J 47r 

The equation above reproduces the well-known kinetic 
equation for impurity scattering 
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(17) 



In this work we will restrict to the case, where en- 
ergies and momenta are restricted even more, namely 
er, e'r, qvpT ^ 1- By expanding ( p^ for small energy 
and momentum and taking the angular average, one finds 
the diffusive equation 



In terms of the quasi-classical Green functions, the charge 
and current read 



/,(x, t) - 2eNo ( J / (P, x) - e(/)(x, t) ] (21) 



2 J 47r' 



j(x,i) = cttA^o / — WFP5tt (P,x). 



(22) 



It is useful to consider the current density in the pres- 
ence of an electric field E(x) — — V(/)(x). By replacing 
in (20) with G^Yi^G^ we express the current density 

as 



j(q,o.)=iey-y ^-G 
(q)G^ ( 



from which we obtain 



uo a 
^2'P+2 

uj q 

e , p 

2'^ 2 



j(x,t)--e^i?iVoV j j5#(P,x) 

= -DWp{x, t) + 2e^DNQE{x, t). 



(23) 



(24) 
(25) 



Within the adopted approximations , i.e., a constant den- 
sity of states and a uniform diffusion coefficient, one ob- 
serves that the current is a linear function of the electric 
field as long as the charge density /9(x) stays uniform. 

We close this section by commenting on Eq. (^l|) . A 
scalar field (/>(x, t) shifts the entire Fermi surface, i.e., it 
affects the Green function at all energies. This is lost in 
the naive substitution of Eq.(^), when the ^-integration 
is extended to ±oo. The second term in eq.(|2l|) is ob- 
tained by taking into account these high energy terms 
correctly. The equilibrium response to a static field is, 
for instance, fully given by this second contribution. For 
this reason the second term is often referred to as the 
"static contribution" to the response, whereas the first 
term is referred to as the "dynamic contribution" . 
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III. QUANTUM CORRECTION TO THE 
CURRENT 



Quantum interference gives rise to corrections to the 
semi-classical expression of the electrical conductivity 
of a metal. The so-called quantum corrections to the 
average conductivity are the weak localization correc- 
tion (WL), the interaction correction in the particle- 
hole channel (EEI) , and the interaction correction in the 
Cooper channel (EEIC). In this paper we will concen- 
trate on the interaction correction in the particle-hole 
channel. For non-linear effects in WL we refer to the lit- 
erature [^,^ . Interactions in the Cooper channel will not 
be considered. This is justified for non-super-conducting 
materials since in that situation the relevant interaction 
parameter scales downwards under the renormalization 
group. 



p + q 



p' + q 



Gi 



FIG. 2. Graphical definition of the diffuson (particle-hole 
ladder). and are the retarded and advanced Green 
functions, which in the general case can depend on external 
electromagnetic fields. 

For completeness we give the expression in the pres- 
ence of external electromagnetic fields. In this case it is 
convenient to go a real space representation where the 
diffuson is defined by the equation 



(30) 



A. Ladder diagrams 



Before calculating the quantum corrections we intro- 
duce the ladder diagrams of repeated impurity scattering 
which will appear at various places in the diagrammatic 
approach. Technically speaking these ladder diagrams 
appear when averaging a product of a retarded and an 
advanced Green function. Here we briefly recall how to 
derive the expressions for the ladder in the absence of 
external fields and without spin effects. The inclusion of 
external fields and spin structure is straightforward and 
one may refers to the reviews on the subject like Ref. [^. 

The diffuson D{q,uj) or particle- hole ladder is found 
by summing the sequence of diagrams shown in Fig.^j: 



RA 



in 



RA\2 



RA 



(26) 



with 



^RA 



1 



2ttNot J (27r)3 
1 - T(-icj + Dq^ 



G^(e + cj,p- 



-q)G^(e,p) (27) 
(28) 



where we have used the condition that lot <ti 1 and 
vfQt <C 1 so that the diffuson reads 



1 



D(q,c^) = - 



1 



(29) 



with Ad — A(x, t -I- rj/2) — A(x, i — r]/2) and = 
(/)(x, t + r]/2) — 0(x, t — rj/2). In these equations t is the 
center-of-mass time, and rj is the relative time, and are 
defined in Fig.^. 



x,t- ?7/2 



x',t'-r?/2 



G' 



FIG. 3. The diffuson in the space/time domain. 

Notice that the external field drops from the equation 
for the diffuson when the relative time rj equals zero. 



B. Interaction correction to the current 

We are now ready to allow for electron-electron inter- 
actions. Interactions will enter the kinetic equation and 
determine the form of the distribution function. We as- 
sume that the distribution function has been determined 
self-consistently via the kinetic equation with the inclu- 
sion of the interaction. We will concentrate then on the 
calculation of the interaction corrections to the current 
density. To do so we need the expression for the Keldysh 
Green function in the presence of interactions. Following 
Refs. [E3 22| we start with the self-energy 



(31) 
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where E™p is the previously defined impurity self-energy 
and 



Sj^(a;,a;') = i ^ J dx2dx3dxidx5Tii,{x5; x, X3) 

i' j' kk' 

X V''''\x5,Xi)Gi,j>{x3,X2)f)',j{x4,;x2,x'). (32) 
The vertex functions are given by 



Among the many contributions to 5Ti we start with the 
Keldysh component of (5I]""p. We denote the correspond- 
ing correction to the current as (5ja, which we determine 

as 



xi,X2) = lij + ^^^^^ Yl j da;id4G'jf (xi, x'l) 
i'j' 

X Ti,ji{x;xi,X2)Gj'j{x2,X2)- 



An analogous equation is valid for F. We recall that in 
the Keldysh triangular representation the "absorption" 
and "emission" vertices differ. A diagrammatic repre- 
sentation of both the self-energy and vertex equations is 
shown in Figs.^ and |5[ 



FIG. 4. The self-energy containing both disorder and in- 
teraction. 



+ 




FIG. 
lines. 



5. Dressing of the interaction vertex with impurity 



The indices denote matrix indices in Keldysh 

space. The bare vertices 7, 7 are local in space and time. 
The structure in Keldysh space is jIj = jfj = 6ij/\/2 
and -ffj = 7j?j = <jfj/V2. From S™p -t- one may de- 
rive a kinetic equation for the system with disorder and 
interaction. The general expressions were already given 
in the seminal paper, Ref. [ p3| . Explicit expressions for 
the various components of F and F in terms of integrals 
like 7]^^ are also given in the appendix of Ref. ||2^ . Fully 
evaluating the expressions in this or that limit remains 
still to be done. Even in thermal equilibrium we are not 
aware of any full self-consistent calculation. 

In the following we will take into account only the non- 
interacting self-energy self-consistently, and restrict our- 
selves to the perturbation theory for the interacting part 
of the problem. The change in due to the interaction 
may be written as 

SG^ = GHj:^G^ + G^ST.'^G^ + G'^ST.^G'^, (34) 

where 6T1 is a sum of the interaction self-energy plus the 
interaction-induced change in the impurity self-energy 



(5S = (5I]""P E 



V 



(5ja(x,t) = -ieD27riVoTV(5S]""P'-^^'. 



(36) 



(33) 



Apparently 5]a is related to the correction to the charge 
density and may be written as 

5j„(x, t) - -e^DN^V I ^<55« (P, x) (37) 



= -DV5p{yi). 



(38) 



Some typical diagrams contributing to 5}a are shown in 
Fig.|. 





FIG. 6. Interaction correction to the current; diagrams 
of the type shown here correspond to 5ja and may be related 
to the gradient of the density. 



The calculation of the correction to the current from 
all the other contributions, which we will denote by Jjf,, 
simplifies due to the following observation: the majority 
of the components of the rcnormalizcd vertices F, F are 
strongly enhanced over the bare value due to the pres- 
ence of a diffusive type of vertex correction. The explicit 
calculation [E2| shows, however, that the diagonal parts 
of F},- and FTTare not renormalized. In the leading order 
of a gradient expansion we can therefore neglect all the 
terms involving these vertices. As a result the expression 
for the various self-energy components are given by the 
combinations 1221 



and 



2j . i i2'~r22'- 21 ' 

TX'^ : Tl,Giiri2V^ 

TX'^ : T\2G22T\2V^ + Tl^Giifl2V^ 



^j.imp,A . G^Ty^^G^/{2TTNoT). 



(39) 
(40) 
(41) 



(42) 
(43) 



(35) 



A diagrammatic representation of the correction to the 
current from YX^^ and YX^^ is shown in Fig.^ 
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FIG. 7. Interaction correction to the current; these dia- 
grams contribute to Sjt- 

Going through the algebra one convinces one-self that 
the contributions from E^'^ and JE""p may be combined 

as 



(44) 
(45) 



where E^^^ and E 



V,K 



refer to the two terms entering 
the expression for E^'^. These combinations of terms 
are taken into account in the diagrams of Fig.^ by com- 
pleting the "Hikami box" as shown in Fig.H. 




FIG. 8. The three diagrams constituting the Hikami box; 
within the here-applied formalism the diagrams are gener- 
ated from the self-energies + E . 

We can now proceed to the explicit calculation of the 
correction to the current. Let us start with the evalua- 
tion of the Hikami box. The first of the diagrams of the 
Hikami box is shown in more detail in Fig.^ 



;Q + q 




q/2 ; t - 



p + qi - q/2 ; t - ?? 



9tt-Mi) 



FIG. 9. Labeling of momenta and times in the Hikami 
box; within the here-applied approximations the retarded 
and advanced Green functions are local in time. 



The evaluation of the diagram amounts to the inte- 
gration over the "fast" momentum entering the electron 
Green function and is accomplished by 



d^p p 
(27r)3 m 



G^(p-Hq/2)G^(p-q/2) 



xG«(p - q/2 + qi)G^(p - q/2 + qi - Q) 



(46) 



The integration is performed under the assumption that 
the three momenta q, qi, and Q are small compared to 
\/l and that the energy in all the four Green functions 
is small compared to l/r. After expanding the Green 
functions to first order in q, qi, Q, the integration over 
p gives (— 67riA*'oT^)Z?(Q + q). Evaluating the other two 
diagrams of Fig.^ analogously and summing up the three 
of them one arrives at {—4,tt\Not^)DQ,. 

After completing the evaluation of the Hikami box, we 
now address our attention to the vertex functions F and 
F. From Ref. pi] we borrow the relevant expressions as 



FL = i=(l-^«^)-i(^«^-K^^^) (47) 
fL = ^(l-0-\ (48) 



so that F22 is simply a diffusion operator, as it is also 
manifest in Fig.^. The explicit space and time depen- 
dence is determined as 

f\^{x]Xi,X2) = --^i:'J't7°(x,Xi)(5(xi - X2)5(tl - 12). 



(49) 



The vertex Fjj has a more complicated structure due to 
the presence of a Keldysh Green function in the kernel 
of the corresponding integral equation. Its detailed dia- 
grammatic representation is shown in Fig.O. 




q + qi 

FIG. 10. The interaction vertex; compare Fig.^. 
The integral rj^^ is given by 



RK 



X G^'^'(p + q/2-Q;p + qi -q/2-Q). (50) 
We replace G^ by G^(-i/T)i^G'^, with 
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(51) 



and integrate over p in Eq. (|o|) to get r]^^^ (x) = Fu' (x) . 
It can be easily shown that rj^^ does not contribute to 
the current; the (relevant part of the) vertex is then found 
as 

r}2(a;;xi,X2) = d?7^I'i-r,/2,i-r,/2(^i' 

X Ft,t_^(x)<5(xi - X2)5(t2 -ti+ v)- (52) 

The last ingredient entering the expression of the in- 
teracting self-energy is the electron-electron interaction 
propagator V^'^. At the level of the approximation we 
are working, it is sufRcient to confine to the standard 
random phase approximation (RPA) . The retarded RPA 
screened Coulomb interaction reads : 

- [ dxidX2V"ix,Xi)U^ixi,X2)V'^ix2,x'), 



the center-of-mass and relative coordinates, F(ei,e2) 
F(gj+g2)/2(ei — £2)- In order to evaluate the frequency in- 
tegrals we assume that the time dependent distribution 
function F^{t) deviates from the equilibrium distribution 
function only at small energies e, i.e. F^(t) — > ±1 for 
e ±00. We then find 



n? 



d 



^(xi,X2)=2iVoT^i?°t,(xi,X2). 



(57) 



Since here the diffuson enters with the relative time 
?7 = 0, the density correlation function does not depend 
on the external fields. Adding the static and dynamic 
parts together and making use of the differential equa- 
tion for the diffuson, one finds 



U''ixi,X2) = -2NoTDdlDlt^{^,,^2). 



(58) 



In the case 
dard expression 



of a uniform system this reduces to the stan- 



(53) 



n^(q,c^) = 2A^o- 



(59) 



where 



V''{x,x')=S{t-t')ey\x~x'\ 



(54) 



from which the dynamically screened Coulomb interac- 
tion is determined as 



and II^{xi,X2) is the retarded component of the density 
correlation function. We write the density correlation 
function as the sum of a "static" and a "dynamic" part, 
-f n''. The "static" part is given by 



V"iq) 



l + VO{q)2NoDqy{~ 
1 -iu) + Dq^ 



Dq^) 



2N, 







Dg2 



(60) 
(61) 



(xi, X2) = 27Vo<5(xi - X2)(5(ti - t2) 



(55) 



as it is seen directly from the quasi-classical expression 
for the electron density in Eq.(^l]). 



62 + (W2 +W3)/2 

4-(^'2+^'3)/2 



£2 - i^^2 +li'3)/2 



FIG. 11. The dynamical part of the density correlation 
function. 

The "dynamic" part is determined from the ladder di- 
agrams as shown in Fig, 
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The expression is given by 



J , f dei dujo 

/ uji uji LO2 yJ2 

^D[e, + -,e,--;e, + -,e2-- 



2 — ia;ii+icj2^' 



UJ2 



UJ2 

~2 



L02 



UJ2 

~2 



(56) 



where we suppressed the spatial indices. -D(. . .) is the 
diffuson and F is the angular average over as de- 
fined in (51). For the frequencies in F we introduce 



The second line is valid, when < V^{q)2N[jDq'^ . In 
three dimensions, where V^{q) — Ane^ /q^ ^ this condition 
reads \oj\ < Dk^. Since the inverse screening length k is 
in a metal typically of the order of the Fermi wavelength 
the approximation is well justified. In lower dimensions 
a more careful analysis is sometimes necessary, see for 
example Ref. 

If we now collect all the pieces of our analysis we may 
come back to the quantum correction to the current. For 
convenience we switch from the energy/momentum do- 
main to the time/space domain. Remember that we have 
neglected the energy dependence of the Green function 
and in the calculation of the Hikami box. This is 
equivalent to approximate these Green functions as local 
in time. The resulting time dependencies are shown in 
Figs.^ and The correction to the current Sj = Sja+Sjb 
is finally found as 

6iai^,t) = -DVSp{x,t) (62) 
(5jb(x, t) = e2nDNoT^ J d?7da;idx2 

; i^t_^,t (x)L>;'_^/2^ ^1 ,«i (^1 ) 

;V5,(xi,X2)(~iVx)i??,,,_,,(x2,x) +C.C.. (63) 



Xj 

X^ 



This current formula is one of the central results of our 
paper. The nice feature is that is valid for arbitrary 
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form of the distribution function and diffuson propagator. 
This will allow us to examine the current in different ex- 
perimental and geometrical setups as well as more general 
questions concerning its physical interpretations. Specific 
applications are discussed in Sections [V ^. The inter- 
acting disordered electron problem is often formulated in 
terms of the field-theoretic non-linear sigma model and 
one may wonders how the presented diagrammatic ap- 
proach is related to it. To this end we explicitly show in 
the appendix ^ at the end of the paper that the same 
expression for the current can be obtained from the field 
theoretic approach of Ref. [pSl . 



IV. GAUGE INVARIANCE 

As a first application of the formula for the current 
we discuss in this section the issue of its physical inter- 
pretation and of gauge invariance. To begin with, we 
believe useful to make contact with our earlier work in 
Refs. [0Jll, 

In Refs. [ prlJT^ the quantum correction to the cur- 
rent was derived working in a vector gauge, A = 
—tE, (/) = 0, under the assumption that the elec- 
tron distribution function had the equilibrium form, 
F(e,x.) — tanh(e/2T), or, equivalently, in the time do- 
main Ftt'{x.) — —\T I sinh[7rT(t — t')], and that we dealt 
with a uniform system with a homogenous charge density. 
In this situation one has immediately that (5ja = and by 



Fourier transforming Eq.(63) with respect to the spatial 
variables, the current formula of [^7|jl^] is reproduced. 

In lljll we showed that in the limit of weak elec- 
tric field our theory reproduces the well-known Altshuler- 
Aronov corrections to the conductivity. At larger fields 
non-linear contributions to the current arise as a conse- 
quence of the nonlocal character of the current formula. 

By working in the vector gauge, the electric field enters 
the equation for the diffuson, Eq.(^0[), via the minimal 
substitution of the vector potential. Due to the quasi- 
classical nature of the equation governing the diffuson, 
one may interpret the non-linear conductivity in terms 
of phases. Let us recall the argument of Ref . The in- 
teraction correction to the conductivity is related to the 
propagation of a particle and a hole along closed paths. 
Pictorially one may think of this as one particle going 
around a closed path, starting for example at i = and 
arriving at t = rj. This particle is also interacting with 
a background particle which is retracing backwards-in- 
time the same closed path. Since the point of interaction 
x(ii) can be anywhere along the path, the particles tra- 
verse the loop at different times. In the presence of a 
vector potential the accumulated phase difference of the 
two paths is ipi — (p2 — e dt'xi • A — e JJ' dt'x2 ■ A. 
This can be simplified using that xi(i) = X2(t) for 
< t < ti and xi(i — -q) = X2{t) for ti < t < ?] 



leading to fi - ip2 = e J^^_^dt'-ki ■ [A{t') - A{t' + rf)]. 
For the particular case of a static electric field described 
by A = — Et, the above given phase difference becomes 
fi — 'f2— e?y(x2 — Xi) • E. This suggests that the inter- 
action correction should be sensitive to a static electric 
field, leading to a non-linear conductivity. 

One may object against this interpretation by observ- 
ing that the vector potential can be gauged away in such 
a way that the static electric field is described by a static 
scalar potential E{y) = — V0(r). A static scalar poten- 
tial, again according to Eq.(^0|), does no longer affect the 
diffuson propagator, so that the argument of the phase 
difference along the two paths cannot be used. Of course 
this does not imply that the current formula is incorrect. 

In fact we may demonstrate explicitly the gauge in- 
variance of the current formula. First one notices that 
(5ja = —D\/p is gauge invariant. For Sjb an explicit check 
is necessary. Given the gauge transformation 



A 



dtX 



(64) 
(65) 



the diffuson and the distribution function transform ac- 
cording to 

Ftr{x) -> F„,(x)exp{-ie[x(x,t) - x(x,<')]} (66) 
Dl,{x,x') ^ D'^^,{x,x') (67) 

X cxp{-ie[x(x, t + |) - X(x, t - |)]} 
xexp{ie[x(x',i'-f |)-x(x',i'-|)]}. 



By applying the above transformation to Eq.(|63|), one 
easily verifies that the function x(x, t) drops, so that the 
expression is manifestly gauge invariant. For the spe- 
cial example of a static electric field with A — — Ei we 
choose x(x, t) — tE ■ x. After the gauge transforma- 
tion the electric field appears in the distribution function 
tanh(e/2r) tanh[(e - n^)/2T], /ix = eE • x but not in 
the diffuson. 

We conclude that although the correction to the cur- 
rent as derived in this paper is gauge invariant, the inter- 
pretation of the non-linear effects depends on the actual 
choice of the potentials A or 0. With E = — 9t A, = 
we would interpret the non-linear conductivity as due to 
dephasing. With E = — V0, A = the reason of the 
non-linear conductivity is attributed to the different lo- 
cal chemical potential felt by the particle and hole. 



V. MESOSCOPIC WIRE 

In this section we use the current formula to analyze 
the non-linear electrical transport in a thin wire. We 
assume that the temperature is low enough so that the 
Drude conductivity is dominated by the impurity scatter- 
ing and is therefore temperature independent. It is well 
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known that in one dimension the electron-electron inter- 
action leads to a l/^/T correction to the conductivity. 
The interesting question to ask concerns what happens 
at larger voltages and which are the relevant length and 
energy scales in the problem. 

For the calculation we need the diffuson propaga- 
tor and the distribution function in the wire. At the 
boundary with the vacuum or an insulator the deriva- 
tive of the diffuson normal to the boundary vanishes, 
(n • V)£'(x, x') = 0. In the case of an infinitely long wire 
with cross section S the solution of the diffusion equation 
reads 



D{x,t) = 



1 1 



f 



T S ^/inDt 



exp[-x^/{4Dt)]. 



(68) 



In the above result we have averaged the diffuson over 
the cross section. In a wire of finite length we impose the 
open boundary conditions along the x axis 



Dtt'{x,x') 



I x,x'—O.L 



0, 



(69) 



corresponding to the fact that an electron arriving at the 
boundary escapes into the leads, where dissipation takes 
place. Therefore it no longer contributes to the phase 
coherent process of quantum interference. The diffuson 
in the finite system is related to the propagator for an 
infinite system according to 

oc 

Dtt'{x,x')^ [D{x- x' + 2nL,t-t') 

n— — OO 

-D{x + x' + 2nL,t-t')]. (70) 

For the actual calculations it is convenient to consider 
also the product of the retarded interaction and the dif- 
fuson, 

iV''D)tt'{^,^') = J dxidiiKf,(x,xi)i?o,,(xi,x'). 

(71) 

For the case of long range interaction this product solves 
the equation 

[-DVlW^DU'M = -^5(x-x')<5(t-t'), (72) 

as it may be seen by comparing with Eq.(|6l|). For a 
one dimensional wire with open boundary conditions wc 
obtain then 



{V'''D)u'{x,x') 



{L - x')x 



L 
1 



{x ~ x')e{x ~ x') 



2DNoT 



(73) 



Besides the boundary conditions for the diffuson and the 
interaction we need the boundary condition for the distri- 
bution function. We assume that the leads of the wire are 



in thermal equilibrium, so that the distribution function 
is given by 



F(e,x) 



=0,L 



tanh 



e ± eV/2 
2T 



(74) 



The distribution function inside a mcsoscopic wire has 
been investigated both experimentally |l9| and theoreti- 
cally 1 26 2^. In the theoretical analysis, in particular, a 
solution of the Boltzmann equation in the presence of dis- 
order, electron-electron interaction and electron-phonon 
interaction has been given. Here we borrow the approxi- 
mate solutions for the the distribution function found in 
Refs. [|6||27|. 



The form of the distribution function depends on the 
various relaxation mechanisms governing the collision in- 
tegral. In the following we first consider a long wire, 
L ^ iph, and then subsequently reduce the length to 
Lph ^ L ^ Lin and Lin ^ L ^ Lt- Here we indicate 
with Lph, Lin, and Lt the electron-phonon, the inelastic 
and the thermal scattering lengths. 



A. Long wire 

In the case of a long wire, L Lph ^ Ljn, the elec- 
trons which traverse the wire scatter many times inelas- 
tically and exchange energy with the environment, for 
example with the phonons. As a result the distribution 
function acquires the equilibrium form with a local chem- 
ical potential and temperature. Our ansatz for the dis- 
tribution function is 



F{e, x) = tanh 



e + eV{L-2x)/2L 
2Te(x) 



(75) 



We assume that the temperature is constant in the bulk 
of the wire, and we also neglect the region near the leads 
where the electron temperature increases from the value 
in the leads to the one in the bulk. The electron temper- 
ature in the bulk may be estimated with standard energy 
balance arguments p5|jT^ ] . For a stationary temperature 
Te, the Joule heating power Pin — cE^ equals the power 
which is transferred into the phonon system, Pout- For 
weak heating, one has Pout = cyAT/rph, where cy is 
the electron specific heat. One obtains then that the 
difference of electron and phonon temperature may be 
estimated as 



AT; 



r-D(e\^/L)Vph/r. 



(76) 



For strong heating, on the other hand, the effective elec- 
tron temperature is of the order of the voltage drop over 
a phonon length 



(77) 
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We are now ready to evaluate the quantum correction 
to the current in the wire as 



dxf{x) 



(78) 



where is the component of the current parallel to the 
wire. Recall that we separated the correction to the cur- 
rent density into two contributions Sj = Sja + <5j&- The 
first of the two terms does not contribute to the correc- 
tion to the current since 



dxSja cx / dx-^Sg'^ix) 
^0 ox 

= Jg^(L)- 5.9^(0) 



(79) 
(80) 



and the boundary conditions impose that Sg^ vanishes 
on the leads. The correction to the current then reads 



/ = 



27rer 
X D 



J dxdxi y"d77Re |Ff_^,t(a;)i^t„,,,t_2,7(2;i) 

?~r,/2.t-3r,/2(2^' a;i)(-i) [6(0:1 (81) 



It is useful to introduce the center-of-mass and relative 
coordinate R — (x + xi)/2, r = xi — x. The last 
term in Eq.(|8l|) above becomes then Q{r) — R/L — r/2L. 
The term R/L vanishes upon integration due to the 
antisymmetry of the r-integral. The assumption that 
the thermal length is much shorter than the system 
size allows furthermore to approximate 8(r) — r/2L by 
6(r), to neglect the boundary effect on the diffuson, 
i.e -C^_r,/2 t-3r;/2(^' ^1) ~ ^(''':'7)j Eind to extend the r- 
integration to infinity. By inserting the distribution func- 
tion we arrive at 



51 = — 27rer 



d?7 



Jo \sinh(7rTe77)^ 
X D{r,ri)sin{eVrT]/L). (82) 



This is equivalent to what is obtained in Ref. |17|. In the 
limit of low voltage the current becomes 



SliV) 



V -4.92 + 0.21 



D{eV/LY 
T3 



(83) 



where 4.92 and 0.21 are approximate numerical fac- 
tors. In order to obtain the full current, one has to 
add the contribution of the Drude leading term, i.e., 
/ = 2e'^DNoSV/L + SI. The low voltage expansion 
applies when the voltage drop over a thermal length is 
smaller than the temperature, eVLx/L < Tg. Since we 
have assumed that Lp^ ^ Lt, the electron temperature 
as a function of voltage rises so fast that the condition 
always holds. Finally we compare the heating and non- 
heating contribution to the non-linear conductivity at low 



voltage. By taking the linear conductivity and the in- 
crease in temperature due to low voltages from Eq. ([76|) , 
we find the heating contribution to the cubic term in the 
current voltage characteristics in the form 



6Ii 



hcatii 



4.92 



e'Lr,, 3 D{eV/Ly 



-V 



7r2 L 2^2 r2(l/Tph) 



(84) 



This has to be compared with the corresponding non- 
heating cubic contribution 



non— heating 



0.21 — 



D{eV/Ly 

r3 



(85) 



One observes that the heating contribution is by a factor 
of the order of Trph larger than the non-heating contri- 
bution. 



B. Intermediate length 

Now we consider a wire of intermediate length where 
L-ph ^ L ^ Lin. One still expects to be near local equi- 
librium although the main mechanism which carries the 
energy out of the wire is not due to the phonons, but to 
the heat flow out of the wire. Under these conditions, 
a temperature profile over the wire develops. The local 
temperature satisfies the equation ES] 



6 {eVy 



which is solved by 



x). 



(86) 



(87) 



The correction to the current at a point x probes the 
wire in a region ~ Lt around x. As a consequence, 
non-linearities in the I — V characteristics arise because 
(1) the temperature depends on the voltage ("heating"), 
and of (2) the non-local character of the current for- 
mula ( "non- heating" ) . It turns out that the first effect 
dominates, whereas the non-heating contribution to the 
non-linear conductivity is only a small perturbation. For 
illustration, notice that the heating becomes strong at 
eV ~ kT. Non-heating non-linearities, on the other 
hand, arise on the scale eVLx/L ^ kT. Since we as- 
sume that L ^ Lt, heating is indeed dominant. This is 
also demonstrated in Fig. 
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where we plot the conduc- 
tance 6G = 6I/V as a function of voltage, while vary- 
ing the system size, L/Lt, and the distribution func- 
tion Ftt'{x). For large L/Lt the linear conductance ap- 
proaches 5// F w ~{4.92e^/7r'^){LT/L)V. For the smaller 
system size L/Lt = 5, the linear conductance is sup- 
pressed, due to the chosen boundary conditions. The 
long dashed lines and the full lines correspond to the local 
equilibrium distribution ( [75| ) , and to the non-equilibrium 
distribution discussed below, respectively. The short 
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dashed line {L/Lt = 5), instead, is obtained by taking 
into account only the heating contribution, i.e. by calcu- 
lating from Eq.(^l|) the linear conductivity and averaging 
over the a;-dependent temperature. 



Eg. (pl|) , as it is appropriate when Lt/L is not very large. 
The linear conductance is - of course - the same as in the 
local equilibrium situation, whereas the non-linear effects 
are slightly larger. 



5h 



heating 



L 



Ax6(j{T{x)){V/L). 



One observes in FigjlJ that the "heating" contribution 
reproduces with a good accuracy the much more com- 
plicated full calculation. For the longer system with 
L/Lt — 200, we do not plot the "heating" curve, be- 
cause in this case it is practically indistinguishable from 
the long dashed one. A slightly larger non-linear conduc- 
tivity is found in the non-equilibrium situation as it will 
be discussed in the next section. 




J ' 



10 

eV/kT 

FIG. 12. Interaction correction to the conductance I /V 
for a mesoscopic wire as a function of voltage. I /V is plot- 
ted in units of (e^ IK)Lt I L. The full line corresponds to 
the non-equilibrium distribution function (^^. The long 
dashed line corresponds to the local equilibrium distribution 
function (^Hj) with the x-dependent temperature. The short 
dashed line [L/Lt = 5) is the non-linear conductivity due 
to the heating contribution only, Eq.(B8h. 



C. Short wire 

In a very short wire the inelastic length may exceed 
the system size L. When one neglects the inelastic scat- 
tering, the distribution function inside the wire becomes 
a linear superposition of the distribution functions in the 
leads and reads 



F(e, x) = [(L - x)F{e, 0) + xF{e, L)] /L. 



(89) 



In the limit L 3> Lt the analytic calculation of the cur- 
rent proceeds as in the case of subsection A. In analogy 



to 



we arrive at 



51 



~27rer 



dr 



drj 



T 



/o JQ vsinh(7rr77) 
xD{r, rj) sm{eVr])r/ L. 



(90) 



The numerical results for the current-voltage character- 
istics in the presence of such a distribution function are 
shown in Fig.O. Notice that we integrated numerically 



D. The Spin-triplet channel 



Up to now we have neglected all the spin effects. As 
demonstrated in appendix the current formula can 
easily be generalized in order to include also the so called 
spin triplet channels. The general equation for the cor- 
rection to the current is given in the appendix. With the 
two-step distribution function ( |89| ) and for L 3> Lt the 
equation for the current reads 

= -AttcDNqt'^ ^ / dr / dry / d^i 

1=0 °° 

2 



T \ r 
■ rr . s\n{eVri)- 
smh(7rJ rj) J L 

X i?(r,ii) — (fi?)(r,r;-ti), 
or 



with 



{fD){r.t) 



1-27^ 



exp [-r^(l 



(91) 
(92) 



and 7* are the interaction amplitudes in the spin singlet 
(7° — 1/2) and spin triplet (7^'^''^ = 7*) channels. In the 
low voltage limit this expression reproduces the standard 
result for the linear conductivity. 



5(T 



-4.92- 



L 



l-3(VT^2^-l + 7t)/7t • (93) 



Fig.|l^ depicts the current voltage characteristics for 
different strengths of the triplet scattering amplitude. 
Again we compare the full result (^Tj) with the simple 
heating contribution, that is the average of the linear 
conductivity over the temperature profile given in (p7[). 




FIG. 13. Correction to the conductance as a function of 
voltage. Notice that for a strong negative scattering am- 
plitude in the spin triplet channel, 7t, the quantum cor- 
rection changes sign. The dashed line corresponds to the 
hot electron regime (heating), while the full line to the 
non-equilibrium situation. 
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For strong scattering in the triplet channel the quan- 
tum correction changes sign. As a function of voltage, 
the quantum corrections are suppressed. For the case 
of a distribution function out of equilibrium, the non- 
heating non-linear contributions are stronger than the 
pure heating effects. 



VI. CONCLUSIONS 

We have calculated the interaction correction to the 
electrical current for a disordered metal out of equilib- 
rium. In order to do so we have extended the diagram- 
matic approach of Rcfs. [^7|jl^. We have first demon- 
strated explicitly the gauge invariance of our current for- 
mula for the current density. In particular, we have dis- 
cussed how the physical interpretation of the non-linear 
contribution may be described differently depending on 
the gauge choice made. In the scalar gauge one may 
think as the particle-hole pair lying on a different chemi- 
cal potential. In the vector gauge, the chemical potential 
of the particle and hole are equal. In the latter case an 
argument based on phase differences leads to the conclu- 
sion that the electric field affects the correction to the 
conductivity. 

We have successively discussed in some detail the cor- 
rection for the case of a mesoscopic wire. We have as- 
sumed that the wire is attached to "ideal leads" (infinite 
conductance), which are described by means of bound- 
ary conditions for the distribution function and the dif- 
fusion propagator. We have distinguished three different 
regimes. First we concentrated on a wire near local equi- 
librium, with a constant electron temperature. Besides 
heating, there is a contribution to the non-linear conduc- 
tivity due to the nonlocal nature of the current response. 

For short wires, where the electron-phonon length is 
longer than the system size, we found that the conduc- 
tance scales with voltage over temperature. Such a scal- 
ing behavior has been recently observed in a nanobridge 
|p9| . However, the quantitative shape of the "scaling 
curve" is not universal. We have found two different 
curves in the hot electron regime, where the inelastic 
scattering length is shorter than the system size, and in 
non-equilibrium. 

We would like finally to emphasize that in this paper 
we have concentrated on the quantum correction to the 
current, and we have left aside a derivation of the charge 
density and the kinetic equation including the quantum 
corrections. This latter task implies the evaluation of 
the quantum correction to the distribution function. The 
evaluation of the distribution function and charge density 
are directly connected as it is clear when one expresses 
the distribution function in terms of the quasi-classical 
Green functions as — g^F — Fg^, and recalls the 
relation of with the charge density. The inclusion of 



quantum corrections into the kinetic equation for disor- 
dered electrons have been considered for the weak local- 
ization case without electron interaction in Refs. pl| , ^ , 
and in the presence of interactions in Ref. [Q. In the 
latter work, the non-linear electric field effects have not 
been included. We notice that the quantum correction 
to the charge density out of equilibrium can be derived 
following the procedure described in this paper for the 
current density in section III. This task is however more 
involved for the following reason. Within the leading or- 
der in the gradient expansion, one finds from ST,™^^'^ , in 
analogy to Eq. (p^) , the identity Spa = Sp. The diagrams 
in Fig.|9|, responsible for the contribution to the current 
denominated 6jb, give in the case of the charge density 
Sph ~ 0. This happens because, the evaluation of the 
Hikami box with the density vertex is zero in the lead- 
ing order of the gradient expansion. In order to consider 
the next-to-leading terms, one has to take into account 
higher powers in the inverse of epT and ql. This requires 
the evaluation of also the diagrams with only one vertex, 
r or r, renormalized by the diffusion pole. Therefore the 
full expression of the electron self-energy is considerably 
more complicated and such a calculation will be more 
lengthy than the one we presented here for the current 
density. This task, although worth to be done, is beyond 
the scope of the present paper. 
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APPENDIX A: FIELD THEORETIC APPROACH: 
THE NON-LINEAR SIGMA MODEL 

The field theoretic formulation of the interacting, 
disordered electron system was pioneered by Finkel- 
stein in the 80's [Q. Recently, this formulation has 
been extended to the non-equilibrium case by means of 
the Keldysh technique by Kamenev and Andreev |2^] 
and Chamon, Ludwig and Nayak for the case of 
normal-conducting metals and by Feigcl'man, Larkin and 
Skvortsov | |35t for superconductors. Gutmann and Gefen 
have also used the field theoretic description to cal- 
culate the current and zero frequency shot noise. 

Given the already extensive literature available on the 
subject, we believe that, rather than repeating again the 
derivation of the non-linear sigma model, it is perhaps 
more useful, instead, to show how to obtain our current 
formula (^3|) within the non-linear sigma model. In this 
appendix we will discuss the spinless version of the model. 
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following Ref. , postponing the spin effects to the fol- 
lowing appendix. In the absence of the electron interac- 
tion, the action of the non-linear sigma model is given 
by 



iSo 



ttNo 



DTiid^Q f + 4iTreQ] , 



(Al) 



where the so-called long derivative is defined by 

5xQ = VQ + ie[A,Q]. (A2) 

The field Q = QH, (x) must satisfy the constraints = 1 
and TrQ — 0. The electron interaction is described by 
the following term in the action 



i5i = -i7riVoTr$„7"Q, 



(A3) 



with 7^ = (To and = dx- The fluctuations of the field <!> 
are related to the statically screened Coulomb interaction 
and given by 

-i($,(x,<)<i>,(x',0) = \val^S{^-^')5{t~t'). (A4) 

In the case of long range Coulomb forces one finds 
V = 1/Nq. Since $ couples only linearly to Q it can 
be integrated out: 



so that the current density reads 

j = ^enDNoiid^QQ - Qd^Qf"). (A9) 

At a first glance, this might differ from what is found 
following Ref. where the current density is written 
as 

j - ie^i?iVo(Tr72 (a,QQ - Qd^Q)). (AlO) 

On the other hand, by comparing with the expression 
( |A^ ) for the Green function, it is seen that Eq.([Al0|) 
sums the Keldysh component and the (21)-component of 
the Green function, whereas ( [A9| ) takes only the Keldysh 
component. Since the (21)-component is zero the two 
expressions are equivalent. 



1. Propagators 

The saddle point approximation for the Q-field 

Q^p-fi^O (All) 



/g-i7rAroTr$„7°Q\ 



exp 



dxdtTr(7iQ„(x))Tr(72Qtt(x)) 



reproduces the Drude-Boltzmann theory. The quantum 
corrections are found when considering the fluctuations 
about the saddle point. We parameterize Q according to 



(A5) 



ue-^/2 e^/^u 



(A12) 



Here the trace refers to the Keldysh space. The ap- 
pearance of the product of terms containing 7^ and 7^ 
stems from the Gx structure of the interaction matrix in 
Keldysh space. In order to make contact with the dia- 
grammatic approach, we express the Green functions in 
terms of the Q-fields. First we observe that pa 



where u characterizes the saddle-point distribution func- 
tion 



1 F 



uazU, 



Gt^(x,x') = 



Gq ^ 



2r 



Q + *a7° 



(A6) 



-1 

and W parameterizes the fluctuations, 

w 



(A13) 



(A14) 



where the brackets (. . .)q.$ indicate that one has to av- 
erage over the flelds Q and $. For the sake of simplicity, 
we drop the subscripts Q, $ in the following. By using 
the condition that the fields Q and $ are only slowly 
varying in space and time, one finds a relation between 
the ^-integrated Green function and the Q-matrix. In 
particular, upon taking the s-wave component of both, 
one finds 



By expanding in powers of the W^-field, one finds for the 
non-interacting action, 5*0, up to the quadratic order 



(2) _ ttA^o 
1*0 - — 



^5ft'(p,x) = (Qtt'(x)). 



(A7) 



y dxd<idi2Wtit2 (5ti + + Ddl) Wt^ti 

+ y"dxdti---diiu;t,*,(9xFt,t;^)iDt.t-^(9,,Ft;tj|. (A15) 

The long derivative d^w = Vw -I- ie[A, w] can here also 
be written as c)x = V — ieAtj^(x) -I- ieAt2(x). The (ujlu) 
correlations solve the differential equation 



In a similar way the p-wave part of the ^-integrated Green 
function is related to Q by 



1 



D{d^QQ-Qd^Q), (A8) 



d_ d_ 

' dh dt2 
2 



5{'K~^)5{ti~t-i)5{t2-U). (A16) 
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After introducing the relative times r] = t2 — ti, t]' = 
ti — ^3 and the center-of-mass times, t = (ti + i2)/2, 
t' — {ts + t4)/2, we identify this correlator with the dif- 
fuson, 



(«7t2ti(x)u;t3tjx')) 



2t 



Dl,MS{rj-i^'), (A17) 



as it may be seen by comparing with Eq.(^0|). Notice 
that the relative time ry is conserved during the propaga- 
tion so that f] = rf . The {ww) correlator is the advanced 
counterpart of the diffuson. (ww) = since there is no 
term proportional to ww in the action. Finally the (ww) 
correlator is given by 



(wt^tiWtsti) = -t^Nq / dt[. ..dt'^{wt^tiWt'^t'J{wt'^t'^Wt^u) 



(A18) 



This correlator is zero in equilibrium, when d^F ~ 0. 
Electron int erac tions modify these propagators. By ex- 
panding in ( |A5| ) each Q-field to first order in W, the 
interacting part of the action becomes 

= -iiTTNof^V J dt{idF - Fiu)uiw -w- FwF)tu 

(A19) 

which leads to non-trivial modifications of all the {ww) 
correlators. Only {ujlo) is not modified by the interaction 
and remains zero. We consider now {ww), which in the 
absence of interaction is just the diffuson. We find 



+ 'nrN^VFt^ti ((wt^t^wtgt J$ - {wt^tiWt^t^)'!-) 

= J^5{^-^')5{ti-h)5{t2-U). (A20) 

To understand the meaning of the above equation, let us 
consider first the limit t2 — > ti. The interaction depen- 
dent term does not drop from ( A20| ) since the distribution 
function is singular in this limit, -Ftjti ~ — i/7r(t2 ~ ^i)- 
Upon using this relation we arrive at the following simple 
diffusion equation 



-(l-A^o^)| 



Ddl 



{WuWt3ti)<i 



ttNo 



(A21) 



To make contact with the diagrammatic calculation, we 
identify this correlation function with the product of 
the retarded interaction and diffusion propagators which 
have been defined in Eq.(|7l|) 



2 2 



{wttm,u)i = {V'^D)ttJ{t3 - U). (A22) 



The propagator with four different time arguments is fi- 
nally given by 

{wt2tiWt3ti)<i> = {wt^tiWtsu) -iTrNoV J dxdt'^dt'^l 
^{wt2t,wt'^t'^) Ff^t'^ {{wt'^tiWt:,u)-s> - {wt'^t'^wt^tjf)] ■ (A23) 



2. Correction to the current 

Due to the fiuctuations of Q there are corrections to 
the charge and current density. Here we calculate these 
corrections by taking into account the Gaussian fluctu- 
ations. In the derivation, we parallel the lines of the 
diagrammatic approach. We begin by separating the cor- 
rection to the current in two contributions Sj = <5ja -I- ^jti, 
where 6ja is related to the gradient of the charge den- 
sity. We then proceed by expressing Sjb in terms of the 
fields uj and u. In the third step we will then explicitly 
include the interaction and obtain the current formula. 
By writing Q — Q^p + SQ, the correction to the current 
reads 



-SQd^SQ -■■■}). 



(A24) 



The dots correspond to the terms which appear due to 
{d-x.Q)Q- The term Jq is proportional to the gradient of 
the charge density 

<5j Jx, t) = -e7rDNoV{SQlfix)) (A25) 
= -D\76p{x,t) (A26) 

and we do not calculate explicitly here. The second term, 
Sjb, contains all other contributions. By collecting the 
various summands, one gets 

Sh = 2Fd^{dQ^^) + {6Q'')d^{2F) 

+ {SQ^^d^5Q'^) + {SQ^^d^SQ^^) (A27) 

Now we expand Q in powers of W and we arrive at 

5ib = cttDNo [F{{d^w)w)^ + {wid^w))^F] (A28) 
= 27re£iiVoRe {{wd^w)^F) . (A29) 

In the absence of interactions 6^b vanishes, so that tak- 
ing the interactions into account is essential. Inserting 
{ww)g, from eq.( [A23| ) we find 



Sib = iireDNoT^ 



d77da;ida;2Re j(x) 

^/2(x,Xi)Ft,^f,_,,(xi) 



xy«,^(xi,X2)(-iVx)i??,,,_,(x2,x)} , (A30) 

which is identical to what we obtained in Eq.(|62|). 

Notice that the long derivative in the current formula 
(A2£) reduces to the gradient in the final result. This 
happens since the two time indices in the relevant Q field 
are always equal, {wttWi^ti) <i> oc ^(^3 — ^4). 
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APPENDIX B: SPIN-TRIPLET CHANNEL 

Until now we have neglected the spin effects. These 
arise from the fact that the diffusion is described by the 
particle-hole propagator, which may occurs in four spin 
states depending on the relative spin of the particle and 
hole. In the case of Coulomb long range forces, the most 
singular contribution comes from the singlet channel that 
we have discussed throughout the paper. In addition to 
the singlet channel, there are three triplet channels which 
also contribute to the quantum corrections to the current. 
In this appendix we extend our earlier work of Ref. | [T^ ] 
to non-equilibrium. 

In order to take into account the spin-dependent inter- 
actions, we follow Refs. [^,0 and start from an inter- 
action which is local in space and time. The fermionic 
action is of the type 

Zl\o 

(Bl) 

where T and r2 are dimensionless static scattering am- 
plitudes, and is an operator for a Fermion with spin 
s. The trace includes integration over space and time 
contour, as well as summation over spin. The interaction 
can be written in terms of the charge and spin densities. 



2Nn 



-Tr 



(B2) 



with p = J2s^s'i's and = Ess'^scrls'^^s'- At this 
point it becomes convenient to introduce the interaction 
amplitudes in the charge (singlet) channel, V = F — r2/2 
and in the spin (triplet) channel, F* = — F2/2. Then 
we decouple the interaction with the fields ($,B) for 
the charge and spin. By going through the steps of the 
derivation of the a-model one finds the following modifi- 
cations of the action 



iSo 



ttNo 



[DTrd^Qss'd^Qs's + AiTrZeQss] (B3) 



D 



Dqp = D/Z 



No A^qp = N„Z 



(B7) 
(B8) 
(B9) 



We are now ready to consider the fluctuations. In analogy 
to what has been done in the spinless case, we introduce 
the charge and spin components of the fields w and iD, 



Ws 



3 



1 



^ i—O 



V2' 



(BIO) 
(Bll) 



where a'^^, = 6ss' and al^, for i = 1,2,3 are the usual 
Pauli matrices. In the quadratic fluctuations of the non- 
interacting action, 5*0, the spin and charge terms decou- 
ple, since the structure is of the type ^i(;ss'(- . ■)ws's = 
w'^i- ■ ■)w^ . The coupling term i5'i is to flrst order in 
W given by 



i5i = -i^iVoV2y Tr \ B'^-i'^u 





-w' 



(B12) 

where for brevity we denoted the scalar field ^ hy B'^. 
We assume that the saddle point, i.e. the distribution 
function and therefore the matrix u, does not depend on 
spin. The interaction field, $, is easily integrated out, 
with the result that at the level of the quadratic fluctu- 
ations, the spin and charge contributions are decoupled 
even in the presence of the interaction. 

The correction to the current is finally found as 
3 



Sib = enDNoY,Mw'y^w')^F 



(B13) 



i=0 



i=0 



y d7ydxida;2Re {F^_^ ^{x) 



^^t-„/2,ti-„/2(^'^l)^*i:ti-')(^l) 



^,,(xi,X2)(-iVx)i?j;'^;i°(x2,x)}, (B14) 



It 



where again i = corresponds to the charge and i 



iSi ^ -mNoTi [{'P^Sss' + B„ • a,,07"Q.'s] , (B4) ^ 3 correspond to spin channels. -Dj;7(x, x') is the 



with 



2iVo 



(B5) 
(B6) 



quasi-particle diffusion propagator in the relevant spin 
or charge channel, which obeys the differential equation 
given in eq.(|o|), with the only difference that the diffu- 
sion constant D is replaced by the quasi-particle diffusion 
constant Dqp. 7' is the dynamically screened interaction. 



7'(q,c.)=7ni + ynViV< 



qpJ 



In (B3) we introduced the factor Z which arises in the 
renormalization of the cr-model |3^ . One observes that Z 
can be absorbed in a redefinition of the interaction am- 
plitudes, the quasi-particle diffusion constant, and the 
quasi-particle density of states according to 



= 7 



-iuj + Dqpq 



(B15) 
(B16) 



-i(l - 2Y)uj + Dqpq^ ■ 

For the explicit calculations it is convenient to consider 
the product of the dynamically screened interaction and 
the retarded diffuson. 
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(rz?0«'(x,x')= y'dxi7,V(x,xi)i^;;';/(xi,x') (B17) 



which solves the diffusion equation 



[(1 - 2y )5t - DVl] ifD')u' (x, x') = ^^(x - x')<5(i - t') 

T 

(B18) 
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